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ABSTRACT 

Recent observations indicate a lower Ha to FUV ratio in dwarf galaxies than in brighter systems, 
a trend that could be explained by a truncated and/or steeper IMF in small galaxies. However, at 
low star formation rates (SFRs) , the Ha to FUV ratio can vary due to stochastic sampling even for a 
universal IMF, a hypothesis that has, prior to this work, received limited investigation. Using SLUG, 
a fully stochastic code for synthetic photometry in star clusters and galaxies, we compare the Ha and 
FUV luminosity in a sample of ^ 450 nearby galaxies with models drawn from a universal Kroupa IMF 
and a modified IMF, the integrated galactic initial mass function (IGIMF). Once random sampling and 
time evolution are included, a Kroupa IMF convolved with the cluster mass function reproduces the 
observed Ha distribution at all FUV luminosities, while a truncated IMF as implemented in current 
IGIMF models undcrpredicts the Ha luminosity by more than an order of magnitude at the lowest 
SFRs. We conclude that the observed luminosity is the result of the joint probability distribution 
function of the SFR, cluster mass function, and a universal IMF, consistent with parts of the IGIMF 
theory, but that a truncation in the IMF in clusters is inconsistent with the observations. Future work 
will examine stochastic star formation and its time dependence in detail to study whether random 
sampling can explain other observations that suggest a varying IMF. 

Subject headings: galaxies: dwarf — galaxies: statistics — galaxies: star formation — galaxies: stellar 
content — ultraviolet: galaxies — stars: statistics 



1. INTRODUCTION 

The stellar initial mass function (IMF), an essential 
ingredient for numerous astrophysical problems, is com- 
monly assumed to be invariant with time and with galac- 
tic properties. This ansatz follows from a lack of evi- 
dence for its variation, des pite searches covering a wide 
range of en vironments (e.g. lKroupall2001l : lChabrierll2003l : 
iBastian et al. 2010), more than from theoretical under- 
standing of the processes that regulate the mass distribu- 
tion of stars. However, several recent studies have ques- 
tioned the idea of a universal IMF in both nearby and 
distant galaxies. 

Considering a few examples in the local universe, 
a truncated and/or steeper IMF has been invoked to 
explain an apparent correlation between the masses 
of clusters and the mas ses of their largest members 
(|Weidner fc Kroupal [2001 though see iLamb et"aI1 120101 
for an opposing view), or between galaxy colors and 
Ha eq uivalent widths fHo versten fc Glazebrookl 120081 : 
iGunaw ardhana et al. 2011), as well as between the ra- 
tio of Ha to FU V luminosity (LHn-/ J^FTiv) and the sur- 



face brightness (iMeurer et al 



■ .20091) or the star forma- 

Lee et al r i2009l : IWeidner et aLllMll) 
these claims remain controver- 



tion rates (SFRs; 
in galaxies. However, 

sial, and alternative interpretations have been offered , 



used to argue strongly for a non-unive rsal IMF (e.g. 
IMeurer et all 120091 : IWeidner etHI 120111 ). Recombina- 
tion lines such as Ha ultimately come from the ionizing 
photons produced primarily by very massive stars, while 
FUV luminosity is driven by less m assive stars (~ 50 Mq 
for Ha, - 10 Mg for FUV, using a lSalDetei<[T955l IMFV 
As a result, the ratio LHa/ipuv can be used as a probe 
for the massive end of the IMF, although the interpreta- 
tion of observations is complicat ed by addition al factors 
(e.g. dust) . Recently , both Meur cr et al.l ()2009l hereafter 
M09) and iLee et al] ()2009l hereafter L09) reported that 
^Ha/^FUV varies systematically with gala xy properties, 
which may imply a variation of the IMF. iBoselli et al.l 
(|2009l hereafter B09) have cautioned that dust correc- 
tion or a bursty star formation history (SFH) can pro- 
duce the observed trends even if the underlying IMF is 
universal, but this idea has not been explored in detail 
with models. While they cannot rule out variations in 
the IMF, L09 also argue that the underlying cause for 
the observed trend is not clear and stochastic effects at 
low SFR need to be explored further. 

Following the se sug gestions (see also 

Williams fc McKed [19971 iCervifio fc Valls-Ga"baudl 



without resort i ng to I MF variations ( Boselli et al.ll2009l: 
Corbelli et all [200l iCalzetti et all I2OIOI: iLamb et al.l 
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200,1 iHaas fc Anderal 120101 : lEldridgell2?)Tli r in this work 
we ask whether the observed iHa/^FUV is consistent 
with a random but incomplete sampling of a canonical 
or modified IMF. To test these hypotheses, we combine 
observa tions from B09, L09 and M09 with models f rom 
slugEI (|Fumagalli et al.l 120111: Ida Silva et aU 120 lit) , a 
novel fully-stochastic code for synthetic photometry of 
stellar clusters and galaxies. 

2. THE OBSERVED GALAXY SAMPLE 
^ http:/ /sites. google.com/site/runslug/ 
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Figure 1. Observed Ho and FUV luminosity in a sample of 
nearby galaxies (L09 red circles; B09 gray squares; M09 blue tri- 
angles). Wc also show (in green), ~ 10* SLUG models for a Kroupa 
IMF with and without clusters (top and middle panel, respectively) 
and for the IGIMF (bottom panel). Analytic predictions for the 
Kroupa IMF and IGIMF arc superimposed (purple dashed and or- 
ange triple-dot dashed lines). White crosses mark the mean of the 
simulated distributions, while the cyan crosses (top panel) are for 
a /c = 1 model with Md^min = 500 M0. 

We compile a sample of 457 galaxies from B09, L09 and 
M09 with integrated Ha and FUV luminosities, together 
with corrections for dust extinction and [N II] contami- 
nation. Due to differential absorption, particular care is 
required when performing dust corrections. Here we fol- 
low the procedure adopted by each author in their stud- 
ies, a choice motivated by the fact that the different dust 
extinctions are consistent across the three samples. We 
refer to the original works for additional discussion. 

Similarly, a careful analysis of the selection biases 
should be a prerequisite for comparison with models. Un- 
fortunately, these data-sets have been assembled adopt- 
ing different selection criteria (a nearly complete volume 
limited sample for L09, a subset of galaxies selected from 
the H I mass function for M09 and a sample of galaxies 
with multi-wavelength observations for B09) that are dif- 
ficult to characterize. We will therefore emphasize only 
those results that are believed to be less sensitive to se- 
lection biases. 

3. STOCHASTIC EFFECTS AND SLUG MODELS 



The IMF describes the relative probability with which 
stars at any mass are formed. At high SFRs, the large 
number of stars guarantees a nearly complete sampling of 
the IMF, and the maximum stellar mass nimax that can 
be found in a galaxy approaches the theoretical upper 



end ' 



of the IMF. At low SFRs instead, due to the 



smaller number of stars that are formed, the probabil- 
ity to find massive stars decreases. The IMF is not fully 
sampled, and mmax ^ "mmux,*- In this regime, stochas- 
tic effects become important to describe the observed 
L-Ra/ Lpuy distribution, regardless of the IMF functional 
form. In addition to finite sampling in mass, iHa/^FUV 
is also affected by finite sampling in time since massive 
stars experience evolutionary phases of short duration 
(e.g. WN or WC phases) during which their ionizing lu- 
minosities can vary significantly. 

Since the majority of s tars is believed to f orm in em- 
bedded star clusters (e.g. lLada fc Ladall2003| ). two addi- 
tional effects should be considered. First, the maximum 
stellar mass in a cluster cannot exceed the cluster mass 
Mcci. Second, at low SFRs, a time-averaged continuous 
SFR results in a series of small "bursts" separated in 
time and associated with the formation of a new clus- 
ter. Provided that stars form with a modest spread in 
ages, gaps in the ages of clusters increase the probability 
to observe a galaxy when the most massive stars have 
already left the main sequence. 

Because it treats these types of stochastic sampling 
correctly, SLUG is the ideal tool to test whether incom- 
plete sampling of a canonical IMF can reproduce the de- 
ficiency of Ha observed in low-luminosity galaxies with- 
out resorting to a modified IMF. For any given SFH, 
a fraction fc of the total stellar population is assumed 
to form in clusters. In this case, SLUG randomly draws 
a set of clusters from a cluster mass function (CMF) 

ip{AIcc\) = M~^^ and populates them with stars randomly 



selected from an IMF 0(m) 



over some specified 



mass interval. Currently, all the stars within a cluster 
are assumed to be coeval. The remaining fraction 1 — /c 
of the total stellar population are formed in the field, 
simulated by randomly selecting stars from an IMF. The 
ensemble of clusters and stars is then evolved with time 
and at each time step the luminosity of individual stars 
is combined to quantify the number of ionizing photons 
and Lpuv foi' the simulated galaxy. 

In this work, we consider two classes of model s . The 
"Kroupa model" is based on a universal iKroupal (j2001l ) 
IMF in the mass interval 0.08 - 120 M© and a CMF with 
(3 = 2 between 20 - 10^ Mq. To test the effects of clus- 
tering we run two sets of simulations, one with fc ~ 
(all stars in the field, similar to the Monte Carlo in L09) 
and one with fc = 1 (all stars in clusters). These ex- 
treme cases br acket all possi b le sol utions, although the 
observations of iLada fc Ladal ()2003[ ) suggest that reality 
is likely closer to /c = 1 than /c = 0. The "IGIMF 
model" is instead based o n the integrated galactic ini- 
tial mass function (IGIMF: iKroupa fc Weidneill2003f ). a 
modified IMF in which TOmax is a function of Med that, 
in turn, depends on the galaxy SFR. The IGIMF imple- 
mentati on follows the minimal-1 and minimal- 2 formu- 
lations (jPflamm-Altenburg et al.l [20071 ) . with a IKroupal 
l|200l IMF and a CMF with /3 = 2, minimum mass of 
20 Mq and maximum mass derived from the SFR av- 
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eraged over 10^ year. In both the Kroupa fc — I and 
the IGIMF model, the IMF integrated over a galaxy is 
modulated by the galaxy star formation and the CMF. 
While in the /c = 1 model the relations between these 
quantities emerge from a pure stochastic realization of 
the IMF and CMF, the IGIMF relics on prescriptions for 
these correlations that imply a truncation in the CMF 
and the IMF. 

For both models, we adopt stellar libr aries with the 
Pado va AGB tracks at so la r met allicity (jBressan et al.l 
19931). the iSmith et all (I2002D implementat ion of 
Hillier fc Milled (I1998D and iLeieune et al.l (|1997| ) atmo- 
spheres. and iMaeder fc Mevnetl ()1988[) winds. We run 
10^ simulated galaxies randomly drawn from a Schechter 
function, si milar to the observed UV luminosity function 
(|Wvder et a l. 2005), with SFR 4.5 x lO'^ - 5 M© yr'^. 
We use 10 time steps of 80 Myr per galaxy to improve 
the statistics. Each of them is independent since the 
80 Myr interval ensures that all massive stars that form 
in one time step are gone by the next one. To further 
improve the statistics at higher luminosity (observations 
are biased towards higher Lpuv), we add 2000 galaxies 
uniformly distributed in logarithmic bins of SFR. Due 
to escape or dust absorption before ionization, only a 
fraction /ua of the Lyman continuum photons produces 
ionizations that ultimately yield Ha photons. To account 
for this effect, we correct the Ha luminosity in the mod- 
els by a factor /hq — 0.95 (BOO). Although the exact 
value for /hq is uncertain ( L09), any choice in th e range 
/h„ = 1 - 0.6 (e.g. B09, iHirashita etldl [200l would 
leave our basic results unchanged. We further emphasize 
that reasonable choices of metallicity, stellar libraries, 
time steps, or even IMF (Kroupa versus Salpeter) yield 
very similar results. 

Slug simulations and theoretical predictions are dis- 
played in Figure [H Stochastic effects are evident in the 
distributions, with simulated galaxies scattered in prox- 
imity to the expected luminositi es for a fully sampled 
IMF. Since previous studies (e.g. iHaas fc Andersl[20T(il) 
have highlighted the importance of the minimum mass 
in the CMF (Md^min) for this type of calculations, we 
display in the top panel of Figure [1] the mean values 
for additional simulations that are similar to the /c = 1 
model, but with A/ci.min — 500 M0. 

4. DISCUSSION 

The main result of this analysis is summarized in Fig- 
ure [T] At Lpuv > 10^'' ergs~^ Hz~^, observed galaxies 
lie close to the value of iHa/ipuv expected for both 
a fully-sampled Kroupa IMF and the IGIMF, but at 
fainter luminosity and lower SFR the data deviates from 
the fully-sampled Kroupa IMF curve towards the fully- 
sampled IGIMF curve. This trend has been taken as 
evidence in support of the IGIMF, but both IMFs are 
subject to stochastic effects. When we properly include 
these using SLUG, we see that realizations drawn from 
the IGIMF are completely inconsistent with a signifi- 
cant fraction of the observed galaxies, particularly at 
Lpuv 10^^ — 10-^^erg s~^ Hz~^. Conversely, realiza- 
tions drawn from a Kroupa IMF span a larger range of 
luminosity and overlap with most of the observed sam- 
ple, with clustering responsible for a further increase in 
the luminosity spread (compare the /c = 1 and /c = 
models). 
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Figure 2. Distributions of Ha luminosity in intervals of Lpuv for 
the observed galaxies (solid black histogram) and for SLUG models 
based on the Kroupa IMF with and without clusters (blue dash- 
triple-dotted and green dash-dotted histograms) and for the IGIMF 
(red dashed histogram). Distributions are centered to (Lhq)i the 
logarithmic mean of the observed Ha luminosity in each bin. The 
two downward triangles indicate the mean Lua for a fuUy sampled 
Kroupa IMF (blue) and IGIMF (red). 



The width in the simulated distributions follows from 
the treatment of mmax and clustering. For a universal 
IMF, TOmax can assume any value up to Jrimax,*, regard- 
less of the SFR. At low SFRs, realizations that lack mas- 
sive stars are frequent and skew the distribution to low 
and low Lpuv- At the same time, realizations with 
massive stars are still possible and some models are dis- 
tributed near or even above the theoretical expectation 
for a fully sampled IMF. The narrower scatter found for 
the /c = model emphasizes that stochastic sampling 
of the IMF alone cannot reproduce the entire range of 
observed luminosity (cf. L09). In the /c = 1 model, 
simulated galaxies that lie at the lowest Ljja for any 
given Lpuv have an excess of older and massive clus- 
ters. This is because at low SFR, a massive cluster rep- 
resents a significant event in the galaxy SFH that, on 
average, increases the time interval between the forma- 
tions of clusters and produces an intrinsic level of bursti- 
ness that results in a wi der luminosity distribution (see 
also Ida Silva et al.l 120111 ) . This effect is obviously am- 
plified in models with Afci,min = 500 Mq, as evident 
from the lower Ha/FUV luminosity in these simulations. 
Note that irregula r SFHs are typical of dwarf galaxies 
ijWeisz et al.l 120081 ) and models with bursty star forma- 
tion may reproduce the observed luminosity distribution 
equally well. Conversely, in the IGIMF theory, at lower 
SFRs, mmax *C TOmax,*- Further, only clusters with low 
mass can be drawn. The narrower mass range that is 
accessible translates in a narrower spread in luminosity, 
and none of the models can significantly exceed the ex- 
pected luminosities for a fully sampled IGIMF (see L09). 

Due to incompleteness, statistical comparisons be- 
tween data and simulations are not straightforward. This 
is complicated by the fact that, at any SFR, both Ha and 
FUV luminosities are subject to scatter and a truly inde- 
pendent variable is lacking. Nevertheless, in Figure[2]and 
Table [H we attempt to quantify the agreement between 
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Table 1 

Summary of the Ha statistics in four bins of FUV luminosity from Figure [2] 



Type=- 




Number'' 


Mean logL 


Ua'^ Fully-sampled logLHa'* 


Dispersion" 


Probability' 


N(o-)s 








(erg s^-^ 


) (erg s-i) 










24.0 < logLpuv < 25.0 


Data 




10 


37.1 




0.64 






Kroupa /c = 


1 


241817 


37.0 


37.7 


0.86 


0.5995 


0.5 


Kroupa fc = 





362593 


37.3 


37.7 


0.55 


0.5515 


0.6 


IGIMF 




64253 


35.6 


36.4 


0.72 


0.0034 


2.9 


25.0 < logLpuv < 26.0 


Data 




86 


38.5 




0.49 






Kroupa fc = 


1 


25571 


38.3 


38.7 


0.55 


0.1073 


1.6 


Kroupa fc = 





24705 


38.5 


38.7 


0.30 


0.1205 


1.6 


IGIMF 




9363 


37.7 


38.2 


0.45 


0.0000 


6.9 


26.0 < logLpuv < 27.0 


Data 




153 


39.5 




0.45 






Kroupa fc = 


1 


4248 


39.5 


39.7 


0.40 


0.3524 


0.9 


Kroupa fc = 





4238 


39.6 


39.7 


0.29 


0.0002 


3.7 


IGIMF 




3576 


39.2 


39.6 


0.37 


0.0003 


3.6 


27.0 < logLpuv < 28.0 


Data 




135 


40.5 




0.35 






Kroupa fc = 


1 


3494 


40.6 


40.7 


0.34 


0.4042 


0.8 


Kroupa fc = 





3562 


40.6 


40.7 


0.29 


0.0090 


2.6 


IGIMF 




3296 


40.4 


40.7 


0.33 


0.0365 


2.1 



^ Type of distribution for which the statistics are listed. 

^ Number of observed and simulated galaxies included in each log Lpuv bin. 

Mean of the log Lua distributions. 

logLna for a fully sampled IMF. 
" Standard deviation of the logLna distributions. 

' Kolmogorov-Smirnov probability P associated to the hypothesis that the observed and simulated distributions 
are drawn from the same parent population. Although the absolute values for the listed probabilities arc difficult 
to interpret given the poorly characterized selection biases, the relative differences between the /c = 1 and the 
IGIMF models support quantitatively our conclusion. 
^ Equivalent number of standard deviations A'^ = \/2 erfc~'(P). 



models and observations by comparing the statistics of 
the Ha distributions in intervals of Lpuv- 

Considering the center of the Ha distributions, the 
IGIMF and fc ~ 1 realizations diverge, moving towards 
lower FUV luminosity. At the lowest Lpuv, the two dis- 
tributions are separated by ~ 2 dex, with the data clearly 
favoring the Kroupa IMF. Figure [D demonstrates that a 
truncated IMF and CMF as in the IGIMF produce sys- 
tematically lower Lua than observed. This is particularly 
evident for the interval 25 < log ipuv < 26 that is char- 
acterized by large enough statistics (~ 100 galaxies) and 
likely not affected by severe incompleteness. To quan- 
tify this claim, we perform a Kolmogorov-Smirnov test 
comparing the observed distribution with 10000 random 
sub-samples of the models, extracted to have a size com- 
parable to the data. We report the results in Table [TJ 
We see that the observed data and the Kroupa /c = 1 
models are generally consistent with being drawn from 
the same parent distribution, while the hypothesis that 
the data and the IGIMF models originate from the same 
parent distribution can be ruled out. The /c = model 
also appears to be inconsistent with the data, though not 
by as much as the IGIMF model. 

Finally, in Figure [3] we compare for each individ- 
ual data point the observed iHa,obs with the simu- 
lated Lua^mod, avcragcd over an interval of FUV that 



is the larger of 0.1 dex and the error on the observed 
Lpuv- This difference is then normalized to dtot — 
V ^obs + '^mod where (Jobs is the error on LHa,obs and 
Cmod is the standard deviation of the models. This quan- 
tity is defined such that a distribution centered on zero 
with dispersion of unity indicates perfect agreement be- 
tween models and data (see right panels). 

In Figure [31 the center of the /c = 1 model (top 
panel) is consistent with the observed galaxies within 
2(Ttot for ipuv < 10^^ erg s~^ Hz~^ and within Scrtot 
for brighter luminosity, where tiobs becomes comparable 
to (Tniod- Conversely, due to the narrower scatter, the 
fc — model (central panel) is only partially consis- 
tent with observations, as previously suggested by the 
KS test. Also, simulated galaxies have higher Lho com- 
pared to the observed galaxies. Based on results of two 
independent tests, we favor the /c = 1 model over the 
/c = simulations. The IGIMF model (bottom panel) 
reproduces instead the observed Lhq only at high Lpuv 
and there is a clear systematic offset between the data 
and the models at Lpuv < 10^'' ergs~^ Hz~^. Although 
the usual notion of probability associated to atot does 
not apply since the Ha distributions are not Gaussian, 
we conclude that observations are better described by 
models based on the Kroupa IMF than on the IGIMF, 
particularly at fainter luminosity and when clusters are 
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Figure 3. Difference between Lhq in individual observations and 
tlie mean of tiic simulated galaxies at comparable LpuVi normal- 
ized to the model standard deviation and observational uncertainty. 
Histograms and Gaussian statistics are shown in the right panels. 



10^5 ergs-^Hz-i, 



included. 

Data are very sparse below ipuv 
where the separation between a universal IMF and a 
steeper/truncated IMF is most evident, precluding us 
from concluding that no variation in the IMF occurs at 
very low luminosity. Moreover, we have not explored cor- 
relations with the galaxy surface brightness, suggested by 
M09 to be the physical quantity related to the IMF varia- 
tion. Also, selection biases arc not well characterized and 
a putative population of galaxies with Xhq < 10^^ ergs^^ 
and Lfuv ~ 10^^ - lO^^ erg s'^ Hz^i would call into 
question a universal IMF. Further we assume that all 
stars born in an individual cluster are coeval, and fu- 
ture work will be required to investigate the effects of 
relaxing this assumption. Conversely, the observed ex- 
istence of galaxies with Lhq ^ 10'^* — 10^^ erg s~^ and 
LpuY ^ 10^^ — 10^^ ergs~^Hz~^ poses a direct challenge 
to the truncation in the IMF and CMF as currently im- 
plemented in IGIMF theory. As illustrated in Figure [TJ if 
the IGIMF model is correct, it should be impossible for 
galaxies to occupy this region of luminosity space. The 
conclusion is strengthened by the fact that we have im- 
plemented only a minimal IGIMF. The standard IGIMF 
would predict even less L^ia, exacerbating the discrep- 
ancy with the observations. Uncertainties in the mea- 
surements, particularly in the dust corrections, together 
with the fact that our models arc highly idealized (sim- 
ple SFH, single metallicity, lack of any feedback), may 
explain the existence of galaxies non-overlapping with 
IGIMF models, but the discrepancy wc found is most 
likely larger than the errors at faint Lho- 



5. SUMMARY AND CONCLUSION 

Using SLUG, a fully stochastic code for synthetic pho- 
tometry in star clusters and galaxies, wc have compared 
the Ha to FUV luminosities in a sample of ^ 450 nearby 
galaxies with models from a universal Kroupa IMF and a 
modified IMF, the IGIMF. Our principal findings are: i) 
simulated galaxies based on a Kroupa IMF and stochas- 
tic sampling of stellar masses are consistent with the ob- 
served distribution; ii) only models where stars are 
formed in clusters account for the full scatter in the ob- 
served luminosity; iii) realizations based on a truncated 
IMF as currently implemented in the IGIMF underesti- 
mate the mean Ha luminosity. 

Based on this result, and since other factors not in- 
cluded in our simulations (e.g. dust, escape fraction of 
ionizing radiation, bursty SFHs) can mimic some of the 
features of the observed luminosity distribution (B09, 
L09), we conclude that present observations of the inte- 
grated luminosity in nearby galaxies are consistent with 
a universal IMF and do not demand a truncation at its 
upper end. 

While we show that the current IGIMF implementa- 
tion provides a poor description of available observations, 
our work is consistent with the fundamental idea behind 
the IGIMF, i.e. that the SFR, the CMF, and the IMF 
jointly produce the observed luminosity distribution of 
galaxies. However, our analysis emphasizes that the cor- 
relations between these quantities emerge naturally from 
stochastic sampling and that a further modification to 
the IMF in clusters as proposed in the IGIMF model is 
not needed to account for the integrated luminosities in 
galaxies. Further, our calculation highlights how clusters 
introduce the level of burstiness required to fully account 
for the observed luminosities, owing to the combined ef- 
fects of the cluster age distribution and the short life-time 
of massive stars. Time dependence might be a key ele- 
ment currently missing in models of the integrated IMF 
in galaxies. 
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